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Abstract 



We present an TV = 2-supersymmetric mechanical system whose bosonic sector, with two degrees 
of freedom, stems from the reduction of an SU(2) Yang-Mills theory with the assumption of spatially 
homogeneous field configurations and a particular ansatz imposed on the gauge potentials in the 
dimensional reduction procedure. The Painleve test is adopted to discuss integrability and we focus 
on the role of supersymmetry and parity invariance in two space dimensions for the attainment of 
O l' integrable or chaotic models. Our conclusion is that the relationships among the parameters imposed 

D , by supersymmetry seem to drastically reduce the number of possibilities for integrable interaction 

i-C ■ potentials of the mechanical system under consideration. 

PACS. 11.30.Pb, 11.15.Kc, 05.45.Ac, 02.30.1k 
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1 Introduction 

The study of integrability in classical and quantum field theories has been developed for quite a time, 
actually since the beginning of the eighties, with relevant results that contributed a great deal for the un- 
derstanding of these theories and, moreover, allowed the improvement of non-perturbative techniques £Q- 
[2]. On the other hand, a number of streams of investigation on chaos has been pushed forward, mainly 
considering spatially homogeneous field solutions and by performing calculations in the framework of 
lattice field theory HU-d]- These studies revealed the existence of chaotic solutions in a considerably 
vast class of gauge theories and, more recently, also in the context of superstrings and supermembrane 
theories |T3 | -|T5 ] . 

Up to now, a detailed analysis relating supersymmetry and chaos, in much the same way as chaos 
is studied in field theories, is lacking in the literature. Close to this issue, we should mention a number 
of attempts to discuss stability and chaos in the framework of brane theories, by concentrating on their 
bosonic sector ^J-^H]. Nevertheless, even in this context, one should put more emphasis on the specific 
role of supersymmetry in the determination of stability and chaos. 
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A similar situation is observed in connection with the investigation of integrable supersymmetric 
theories, where the integrable or non-integrable character is ascertained, without however highlight- 
ing the mechanisms or those specific properties of supersymmetry which work in favour, or against, 
integrabilitv[TB ] -PD l . 

Our work sets out to tackle this issue, that we believe should be more manifestly worked out. To 
pursue an investigation focusing on the role of supersymmetry in connection with integrability and chaos, 
we propose to start off from a supersymmetric mechanical system, rather than a field-theoretic model. 
The system we choose to work with is built up as the N = 2-extended supersymmetric version of 
a dimensionally reduced SU(2) Yang-Mills theory that arises when spatially homogeneous fields are 
considered and a particular ansatz on the gauge potentials is adopted in the dimensional reduction 
scheme so that only two degrees of freedom survive[2I] in the mechanical limit. We also devote special 
attention to the role of parity symmetry, since we assume the latter is an invariance of the interactions 
involved in the systems we shall be considering. Our analysis of integrability shall therefore rely on our 
considerations on supersymmetry and parity invariance. They dictate special conditions in the space of 
parameters so that, instead of having to take by decree special choices of these parameters, as it is usually 
done, we invoke these two invariances to naturally restrict and select possibilities in parameter space. As 
a matter of fact, we anticipate that parity may appear in two versions for planar systems, and this point 
shall be suitably taken care of here. 

Our paper is organised as follows. In Section 2, we propose a general 2-dimensional purely bosonic 
model with parity symmetry and we identify the cases of integrability. Next, the N — 2-supersymmetric 
extension of the model is written down in Section 3. The complete bosonic sector, now enlarged by the 
presence of two supersymmetries, is discussed in full details in Section 4, where we pay due attention to 
the role of parity and we pick out Painleve test as a criterium to infer about integrability. In Section 5, 
we reassess the question of the integrability for the bosonic sector of our N — 2-model, but now taking 
into account the constraints dictated by parity whenever it is imposed also to the fermionic interactions. 
A very restrictive class of potentials comes out that fulfills integrability. Finally, in Section 6, we present 
our Final Discussions and we draw our General Conclusions. 

2 The ordinary bosonic model with considerations on parity 
symmetry 

We assume the most general fourth-order polynomial potential for two degrees of freedom described by 
the variables x and y: 

V = dx 4 + C 2 y 4 + C 3 x 3 y + C 4 xy 3 + C 5 x 2 y 2 + C 6 x 3 + C 7 y 3 + C 8 x 2 y + C 9 xy 2 + C 10 x 2 + C n y 2 + C 12 xy. (1) 

It may be considered as a sort of protopotential used to build up a general non-supersymmetric poly- 
nomial potential up to fourth order. We are bound to fourth order because we have in mind mechanical 
models derived from Yang-Mills theories and these, as we know, display self-interaction vertices for three 
and four potentials. Since we are interested in realistic models, we impose parity symmetry which is 
respected by mechanical and electromagnetic models. We shall not be dealing with models coming from 
chiral gauge theories. 

To implement parity in the model, we have to consider that there are two possibilities, since we are 
in a 2-dimensional space: 

x — * — x 

x-parity : , (2) 

y-> v 

or _^ 

2/-parity : . (3) 

y^-y 

In the first case, the resulting potential is 

V = C lX 4 + C 2 y 4 + C 5 x 2 y 2 + C 7 y 3 + C 8 x 2 y + C w x 2 + C llV 2 . (4) 
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This potential looks like the sum of two well-known potentials: 
a quartic potential (Yang-Mills-type) 



V YM = Ax 2 + By 2 + ax 4 
which is known to be integrable in the following cases 22 



2„,2 



+ by + dx y 



(5) 



a) A= B, a = 6, d — 6a. 


that in our case is Cio = Cn,Ci = Ci,C^ = 6C\. 


b) A, B, a = b, d = 2a. 


that in our case is Cio, Cn, C\ = C 2 , C5 = 2C\. 


c) A = 4B, a= 166, d = 12a. 


that in our case is Cio = 4Cn, C\ = I6C2, C5 = 12Ci 


d) A = 4B, a = 86, d = 66. 


that in our case is Ciq = C\i,C\ = 62,65 = 6C\. 


e) d = (trivial) 


that in our case is C5 = 0. 



and the Henon-Heiles potential: 

Vhh 



By 2 



ax y- -by , 



(6) 



that exhibits well-known integrable cases|22| 



a) A = B, a = -6. 


that in our case is Cm = Cn, C7 = ^Cg. 


6) A, B, 6a = -b. 


that in our case is Ciq, Ch, C7 = 2C%. 


c) 16A, B, 16a =-6. 


that in our case is Cio = 16Cn, C-j = ^-Cg. 


d) a — (trivial) 


that in our case is C$ = 



3 The supersymmetric model 

Now, we shall consider an N = 2 supersymmetric mechanical model[2*3^. defined as follows. The two 
(complex) Grassmannian parameters of the superspace will be denoted by 6 and 9. The two real coordi- 
nates of a planar particle, x and y, are the bosonic components of the superfields coordinates, which are 
given by 

(7) 



x(t, e, 9) = x (t) + et 7l A(t) + At(t) 7l e - -e^ 3 &h(t), 



and 



with: 



Y(t,9,9)=y(t) 



e 



A = 



(8) 



(9) 



where all the A's and ^'s are Grassmannian variables. The 7j's are the Dirac matrices corresponding to 
the two-dimensional Euclidean space under consideration and they may be chosen so as to coincide with 
the Pauli matrices: 7,; = <7j and 73 = —47172 = 03. 6 is Majorana spinor, which, in this particular repre- 
sentation of the 7-matrices, takes the form given in @, where the "bar" stands for complex conjugation. 
On the other hand, A and H are Dirac fermions. Therefore, Eqs. J7HSJ yield: 



and 



X 



Y 



(Ai - A 2 ) - 6 (Ai - A 2 ) - 



(10) 



(11) 



It is noteworthy to remark that it is precisely the combination (Ai — A2) the one that carries the 
fermionic degrees of freedom of X. On the other hand, as for Y, its spinorial degress of freedom are all 
located in (£1 — £ 2 ) ■ 

The supersymmetry covariant derivatives are as below: 



D 
D 



de- 
do- 



iOdt 
iSdt, 



(12) 
(13) 
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they satisfy: 

D 2 = (14) 

D 2 = (15) 

{D,D} = -2idt. (16) 

The super-action to be considered contains, besides the kinetic terms, the most general superpotential, 
up to third order in the superfield coordinates (this implies a fourth-order potential in terms of the physical 
coordinates) 



S = 



J dtdddO jy [DXDX + DYDY] +U(X,Y)^, (17) 

where the first term gives rise to the kinetic terms and the superpotential U {X, Y) is assumed to be given 
by: 

U(X, Y) = k x X 2 Y + k 2 XY 2 + k 3 X 2 + k 4 Y 2 + k 5 XY + k[X 3 + k' 2 Y 3 , (18) 

the fc's being arbitrary real constants. Since the term in XY may be canceled out by means of a proper 
linear transformation (a rotation in the X-Y plane), then the constant fc 5 may be set as zero, k§ = 0, 
without loss of generality. Similarly, the terms linear in X or in Y were not considered, since they may 
be eliminated by a translation redefinition, X' = X + const and Y' = Y + const'. The equations of 
motion may be used to eliminate the non-dynamical degrees, of freedom fj, and, thus, the super-action, 
S = J dtL , yields the following Lagrangian where quartic terms in the potential are present: 



Mx ,„/V< — i t+ , ; , ; \ k 2 + 9k[ 2 4 fcf + 9fc 2 2 4 



6k\k[ + 2k\k 2 3 6k 2 k' 2 + 2k\k 2 ,, 

77 x y — xy + 

M y M 

2k\ + 2k\ + 3k[k 2 + Sk^ 2 2 Qk 3 k[ 3 6k 4 k 2 3 

M X V M~ X M~ V + 

4hk 3 + 2fcifc 4 2 4fc 2 fc 4 + 2k 2 k 3 2 2k 2 2 2kj 2 

x v XV -X V + 

M y M U M M J 

+ [2ik x (Ai?j - Aifa - A 2 6 + A 2 6 + Ai6 - Ai6 - A 2 6 + A 2 ? 2 ) + 

- 2k 2 (^6 ~ 66 - 66 + 66) - 6K (AiAi - AiA 2 - A 2 Ai + A 2 A 2 )] z 
+ [2ik 2 (AiCi - Ai^ 2 - A 2 £i + A 2 £ 2 + Ai6 - Ai6 - A 2 6 + A 2 £ 2 ) + 

- 2fci (AiAi - AiA 2 - A 2 Ai + A 2 A 2 ) - 6k' 2 (a6 - 66 - 66 + 66)] V 

- 2k 3 (AiAi - AiA 2 - A 2 Ai + A 2 A 2 ) - 2fc 4 (66 - 66 - 66 + 66) ■ (19) 

In the next sections, the integrability conditions for this Lagrangian will be discussed, and the influence 
of supersymmetry and parity invariance shall be highlighted. 



4 The bosonic sector and its integrability. 

The direct application of the the Painleve test (for a short review, see Appendix [5J directly to the 
bosonic sector is not actually a good procedure, for the resolution of the systems that appear in the 
analysis becomes very complex. 

In this section, we shall take into consideration the observation that the original model is not invariant 
under the two classes of parity transformations. This may set a more formal framework. 

So, in a first attempt, we will impose parity symmetry, that is a discrete symmetry, only to the bosonic 
sector of the theory and after that we shall check how the constraints imposed by this invariance affects 
the integrability of the model. 
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Adopting invariance under x-parity, we have the following constraints on the coefficients of the po- 
tential: 

C 3 = (6fak[ + 2fak 2 ) = 0, (20) 
C 4 = (6k 2 k' 2 + 2fak 2 ) = 0, (21) 
C 6 = (6k 3 k[) = 0, (22) 

C 9 = (4k 2 k 4 + 2k 2 k 3 ) = 0, (23) 

where the C"s above are the coefficients of the bosonic sector of the original potential for which the parity 
symmetry is broken. 

4.1 Parameters surviving the parity constraints 

Solving the system of conditions for fa 1 k 2 ,k' 11 k 2 ,k 3 e k 4 , we obtain as solution the following possibilities: 

{k[ = k[, k' 2 = k' 2 ,fa = k 4 , k 3 =0,k 2 = 0, fa = 0}, (24) 

{fci = fa, k 2 — k 2 , k 4 — fa, k 3 — k 3 , k[ — 0, k 2 = 0}, 

{k' 2 =0,k[ = k[, k 2 = k 2 , k 3 =0,fa = 0, fa = 0}, 

{k[ = k[,fa = fa,k 2 = -3k[,k' 2 = -l/3fa,k 3 = 0,k 4 = 0}, 

{k' 2 = 0,k 2 = k 2 ,k 3 = k 3 ,k[ = 0, k 4 = -l/2k 3 ,fa = 0}. 

To study the consequences of these solutions we shall present in the next subsection the Painleve's 
test (see Appendix 0) that have been very used in the search for integrable systems for being algorithm 
and with wide application. 

4.2 Applying the Painleve test 

For the first case {k[ = k[, k' 2 = k' 2 , fa = fa, k 3 = 0, k 2 = 0, fa =0}, we have the following potential: 

Potl = 9 k£ 4 fa_ 3 + 2 M 2 (25) 

2 M 2 M y M y M y y ' 

Applying the Painleve test, we obtain four branches referring to the uncoupled systems that survive 
the test. 

For the second case {fc 2 = 0, fa = 0,k' 2 = k' 2l k 4 = k 4 ,ki = fa,k 3 = fc 3 }, we have the following 
potential: 

Pnt - 1 M r 4 4- 9 fc 2 „4 _i_ ( 2fc l + 3fclfcg) o 2 , (9f .\ 

P ° t2 ~ 2M X + 2M 9 + M X V + (26) 
k> (4fak 3 + 2fak 4 ) 2 kJ fe| 

with dominant potencies: 

a 1 =-l,o 2 = -l (28) 
and four branches with the following expressions for the resonances: 

_^ ^ (2fa-3k> 2 ) ^ (3k' 2 + fa) ^ 
fa fa 

that will show integer resonances if we set fa, k' 2 = n^fa where n = { — 1, 0, 1, 2}. 
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For the case n = — 1, it is not possible to determine the resonances. 
For the case n = 0, we have the following potential: 

p * -\p+ 4**« 2 + (4fcA M t ' t4> * 2 » + 4* 1 + 4* <*> 

It docs not pass the Painleve test because there appears a compatibility condition that cannot be 
fulfilled: 

-4iV2(18k%xl - 5fc| - 4fc 3 fc 4 ) = 0. (31) 

This equality is indeed satisfied if hi, k% e fc 4 = 0, but this cancels out the potential. 
For the case n = 1, we have the following potential: 

PoU = ^x 4 + + 3§zV + 2fc 4 |V + ( 4fc i fc 3 + 2fc 4 fci) 2 2 fc| 2 2 fc| 2 (32) 

And now, we obtain four branches with the following resonances: 

-1,1,2,4, (33) 

but with the following compatibility condition: 

-2(-k 4 + k 3 )M = 0, (34) 

to be verified in the resonance j = 1 of the first and of the second branch. Setting fc 3 = fc 4 , the system 
becomes compatible and it passes the Painleve test with only two branches and with the potential now 
written like bellow: 

Pot5 = l^L x i + IM/ + A 2 2 + 2 Ml y 3 + 6^*1^ + 2 $x 2 + 2% 2 , (35) 
2M 2M y M M M M M y ' v ' 

with dominant potencies: 

a\ = — 1, a 2 = -1. (36) 
The values of the resonances for the two branches are: 

-1,1,2,4, (37) 

and for the first branch the coefficients of the dominant terms are: 

1 iM 1 M 

X ° = 2lT'^2 Z V (38) 

For the second branch, the coefficients read as follows: 

HM 1 M 

x ° = -2^ y ° = -2 l i;- (39) 

In the first branch, the arbitrary coefficients are: 

Vi, y 2 and (40) 

and the arbitrary coefficients of the second branch are: 

Vn y 2 and x 4 , (41) 

reminding that the variable to is the fourth arbitrary quantity corresponding to the resonance —I. 
So, the system is of fourth order and possesses four arbitrary coefficients; therefore, it is integrable. 
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For the case n = 2, we have the following potential: 
It was not possible to determine the dominant terms. 

For the third case {^3 = 0, k\ = 0, k 2 — k 2 , k\ = 0, k[ = k[, k' 2 = 0}, we have the following potential: 

Potl = gM!^ + 1M 4 + (2fcl + 3fclfc 2 ) 2 2 (43) 

2 M 2 M M tf ' v ; 

the expressions for resonances in this case are: 

(3fcj + fc 2 ) (-2fc2 + 3fcD r . 

' ' fc 2 ' fc 2 1 ' 

that will show integer resonances if we set k 2 , k[ — n^k 2 where n = {—1, 0, 1, 2}. 
For the case n = — 1, the system passes the test with the following potential: 

it 2 ip p 
P * = 5^ 4 4l» 4 + |-V, (45) 

with dominant potencies: 

en = -l,a 2 = -1 (46) 
and the values of the resonances for the two branches: 

0,-1,3,4. (47) 

For the first branch, the coefficients of the dominant terms arc: 



for the second branch, they are: 

X = -V± mi t y = y (49) 

k 2 

In the first branch, the arbitrary coefficients are: 

y ,x 3 and y 4 , (50) 
and the arbitrary coefficients of the second branch are: 

y , x 3 and y 4 . (51) 
For the case n = 0, the system does not pass the test with the following potential: 

because there appears the following compatibility condition: 

\8iklyj = 0. (53) 

This equation is satisfied only if k 2 equal to zero, but in this case the potential vanishes, and this is 
not an interesting situation. 
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For the case n = 1, the system passes the test with the following potential: 

1 k 2 1 k 2 k 2 

™»-5^ + 5^ + »^V. (") 

with the same dominant potencies: 

ai = — 1, a 2 = -1, (55) 

and resonances: 

-1,1,2,4 (56) 

in the two branches. 

For the first branch, the coefficients of the dominant terms are: 

l.M l.M 

(57) 

and, for the second branch, the coefficients appear as bellow: 

1 M 1 M 

X ° = 2 l k- 2 > y ° = 2 l k;- (58) 

In the first branch, the arbitrary coefficients are: 

x 17 x 2 and x 4 (59) 
and the arbitrary coefficients of the second branch are: 

y 1 , y 2 and x 4 . (60) 
For the case n = 2, the system does not pass the test with the following potential: 

because it was not possible to determine the dominant terms. 

For the fourth case {k[ — k[, k' 2 = — l/3fci, fc 2 = — 3^, fci = k\, = 0, k± = 0}, we have the following 
potential (quartic): 

with k[ =k\. The resonances in this case are: 

0,-1,3,4 (63) 

for the two branches. 

For the first branch the coefficients of the dominant terms are: 



1 v/(-10M 2 - 100&?wf) 
X ° = T,~ fa ^o=.o; (64) 



and, for the second branch, they are: 



_ 1 V(-lOM'-lOOfcjfrg) _ 

10 fa (b5) 

In the first branch, the arbitrary coefficients are: 

y , x 3 and y 4 (66) 



8 



and the arbitrary coefficients of the second branch are: 



y , x 3 and y i . 



(67) 



For the fifth case {k' 2 = 0,k 2 = k 2 ,k 3 = k 3 ,k[ = 0, k± = — l/2k 3 ,k\ = 0}, we have the following 
potential: 



This potential does not pass in the Painleve test because the following compatibility condition appears: 



that is only satisfied if fc 2 = fc 3 = 0, and this eliminates our potential. Therefore this case does not pass 
the Painleve test. 

As this potential is of the quartic type, it is easy to verify that the result of this Painleve analysis is 
in agreement with the conditions of integrability for this potential type. 

5 The integrability of the bosonic sector with parity considera- 
tions for the complete model. 

As verified in the previous section, by imposing parity to the bosonic sector, the task of finding integrable 
cases became less arbitrary, in that the choice of the coefficients in the terms of the potential was guided 
by the argument of parity invariance. In spite of that, it was still necessary to fix by hand the values of 
some parameters when applying Painleve test to recover the integrable cases we have listed previously. 

In this section, we shall impose the parity symmetry not only to the bosonic sector but also to the 
fermionic interactions, and we shall verify to which extent the constraints on the parameters are able to 
turn the model integrable without the need of fixing arbitrarily parameters in the Painleve test. 

5.1 Two-component formulation of the fermionic sector 

Since the model is classic and non-relativistic, and defined in a two-dimensional Euclidean space, E 2 , the 
covariance group is SO(2). We adopt the representation below for the Clifford algebra: 
Therefore, we adopt: 



Pot 13 




(68) 



-3ik 2 + I8ik 2 y 2 = 0, 



(69) 




(70) 




(71) 



73 = -i7V = &z, 



(72) 



such that: 



{7\7^}=2<Pn, 



(73) 



{7\7s} =0. 



(74) 



For a general spinor 




(75) 



the action of SO (2) is as given below: 

where u is the rotation angle; therefore is invariant. 



(76) 
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Now, we try to identify x-and-y-parities in the spinorial space. 
To do that, we start off from the Dirac equation: 



i^d^ + i 1 2 d v t> = 0, (77) 

to which we impose ir-parity symmetry: 

*(t;2) P , &(t';x') = (78) 

= RV(t';-x',y'), 

where R represents the parity matrix in the spinor space: 

= -Ry 1 (79) 
7 2 i? = i? 7 2 . 

Then, our parity matrix may be chosen as 

R = 7 2 (80) 

and, thus, 

*'(t'; x') = jH(t;x). (81) 

So, all spinors, up to a phase factor, transform under parity by means of the 72-matrix. 
Considering the other possibility, that is, the y-parity, one can readily check that parity is represented 
by the 71-matrix: 

p{ X ^ X (82) 

t y -> -y 

(83) 



5.2 The integrability with the parity constraints from the fermionic Sector 

To include the constraints dictated by x- or y-parity symmetry for the complete (bosnic + fermionic) 
model, we propose to actually carry out the analysis directly in terms of the superfields J7J) and ©■ 
Rather than following the lengthy procedure of considering all the terms of the component-field action, 
we propose to work without quitting superspace. 

The action of the x-parity on the superfields is given by 

X -> -X and Y -> Y, (84) 

provided that 

7 2 e, 
72A, 

- 72^, 
fu 
- h- 



e 

A 



h 
h 
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With these parity assignments to the fermions and auxiliary fields, the superfield coordinates transform 
under parity exactly as above. Moreover, by virtue of the specific choice of 72, we have that parity acts 
on d6, d9 and the covariant derivatives as below: 

D -> -iD, D -> iD; 
d9 -> idd, d9 -» -id0. 

With all the prescriptions, the volume element dtdOdO picks a minus sign. This means that the kinetic 
terms are naturally invariant, but parity symmetry of the potential sets 

ki = k 3 = fc 4 = k' 2 = 0, (85) 

with ki and k' x non-vanishing. 

These parameters constraints are the same as the third set we found when only the bosonic sector 
was considered and we found only two integrable cases: Potentials 8 and 10 the we rename now as below: 

p «^-'Am* + \m* + m**> (86) 

and 

So, from all integrable cases found when we considered only bosonic sector, only the two potentials 
above preserve x-parity under complete model consideration. 

On the other hand, if we contemplate y-parity symmetry for the whole model, we have that 



X and Y -» —Y, 



provided that 



e -> 7i0, 

A -> - 71 A, 

h —> - h, 
h - h- 

Also, D -iD, D -» iD, dd -> id0 and d9 -» -idS. 

So, as in previous case, y-parity invariance is ensured only for those superfield monomials that change 
sign under parity. This then impose: 

k 2 = k 3 = ki = k[ = 0, (89) 

while k\ and k' 2 are the only coefficients compatible with y-parity invariance. 

These constraints on the parameters correspond to only one set of solutions that is found when only 
the bosonic setor is considered in connection with the y-parity, in a similar way to what happens for 
cc-parity. There are only two integrable cases that we shall present below: 

Po/.,,,,,,.,-,, = T - 7 , ' + '- r : V r , r ,-V. (00) 

and 

Pot susv2 - v = ^^a; 4 + + 3^x z y z . (91) 

So, from all integrable cases found when only the bosonic sector is considered, only the two potentials 
above preserve y-parity if whole model is analysed. 



1 k 2 

-—x 4 - 
2M 


1 k 2 

f -— y A - 

2M y 




1 k 2 

2M 4 


1 k 2 

2M y 


k\ 
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6 Final discussions and general conclusions. 



Along the previous sections, we carried out an integrability analysis of the bosonic sector of the super- 
symmetric model and we verified the appearance of integrable cases for both coupled and non-coupled 
systems. 

The coupled cases turn out to be classified into two types: a quartic potential and a potential that is 
functionally the superposition of a quartic and a Henon-Heiles potential. 

Contrary to the situation where we impose parity symmetry to the complete action (bosonic and 
fermionic interactions) and the generated potentials come out totally integrable, without the need of 
setting integrability constraints, the case in which parity symmetry is imposed only to the bosonic sector 
yields integrable potentials only after we take into account the constraints that appear in the course 
of Painleve analysis. This means that, if these constraints arc not fulfilled, we will be dealing with 
non-integrable potentials and therefore with the possibility of chaos. 

For the cases where the potentials have a quartic form, there is no need to go through a chaos analysis 
for this has already been discussed in the literature we have previously referred to. 

The cases for which the potentials are given by the superposition of a quartic and a Henon-Heiles form 
are under consideration and, in a forthcoming work, we shall report the results of a complete analysis 
|24| . However, in this section, we shall give an example to illustrate how this type of non-integrable 
potential admits order-chaos transition by using the potential of number 6 of Section(4.2): 

Pot6 = \ 1 $-x i + 2§y 4 + S-x 2 y 2 + 4fc 4 |V + ( 4fc i fc 3 + 2fc 4 fci) 2 |[ 2 fcf 2 

21/ M y M M y M M M v ' 

For this purpose, we make use of Lyapunov characteristic exponent (LCE) and phase portraits |25ll2f)| 
and |27j . The Lyapunov exponent is a usefull tool to quantify the divergence or convergence of initial 
nearby trajetories for a dynamical system. In a chaotic system, there is at least one positive Lyapunov 
exponent, defined as 



r , di(t) 
Oi = lim m — — — 

t^oo d% (0) 

where di (t) is a deformation measure of the small hypersphere of initial conditions in the phase space 
along the trajectory. The asymptotic rate of expansion of the largest axis is given by the largest LCE. 
By phase portrait we mean a graph of the dynamical variables in phase space that is used to provide a 
qualitative insight of the dynamical behavior of the system under study. The accuracy of our computation 
was verified by checking if the Hamiltoniam was conserved during the simulation. 

Fixing k\ = 10, M = ks = fct = 1 the potential acquires the following form: 



V := 50x 4 + 200y 4 + 400a:V + + W% V + 2x z + 2y z 

We calculate de largest <Ji and its respectives phase portraits, and we present two cases for the same set 
of parameters fixed above, but with different initial contions. First withpi(O) = 0.1,p2(0) = 0.1,(?i(0) = 
0.1,52(0) = 0.0, Energy=0.035; it presents regular behavior (see Figures 1 e 2 below). 
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The second case is given by Pl (0) = 0.1,p 2 (0) = 0.1, ft (0) = 0.1,ft(0) = 0.18, Energy=0.78; it 
presents chaotic behavior (see Figures 3 e 4 below). 




Figure 3 
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q 2 (0)=0.18 




It is perhaps noteworthy to say that a potential of the form (x 2 + y 2 ) can be generated by considering 
the work of Refs. EHJ-IS], where a classical mechanical Yang-Mills system with four degrees of freedom 
£i is studied, with Hamiltonian given by 

h = \ {p? + p! + Pi + p!) + ^{& + e 2 + £ 3 2 + &) ■ (93) 

Indeed, if we adopt hyperbolic coordinates on the (^15^2) and (£3; £4) planes, 

£l = 7"! COsh#! 

£2 = T\ sinh^! 
£3 = r 2 cosh 6 2 
£,± = Ti sinh 8 2 , 

the Hamiltonian becomes: 

H = \ {Pi + P\ - ]fl ^Pl) +^(rl+ r 2 f . (94) 

9\ and O2 are ignorable coordinates, so that their corresponding momenta are integrals of motion. 
If Pg = Pq — 0, no negative contributions to the kinetic energy are present and we get an effective 

two-degrees-of-freedom system with the potential of the form (a; 2 + y 2 ) 2 . Besides the discussion of 
integrability and classical chaos, studies of quantum chaos using this potential have received a great deal 
of attention in the literature |32|-p?7]. This observation may be a good support in favor of the results we 
have got on the integrable potentials produced after parity symmetry has been imposed for the complete 
supersymmetric model. 
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A Painleve test 

The Painleve test [23 — [2HJ establishes if a system of ODEs exhibits the Painleve property. 

An ODE has the Painleve property if its solutions in the complex plane are single-valued in the 
neighborhood of all its movable singularities. Given a differential system 

Lj(ui, ua) = with i,j = l,...,n , (95) 

we assume a Laurent expansion for the solution 

oo 

Ul {t) = {t-t ) a >Y, u ^{t-to) k , (96) 

k=0 

with 

u i,o 7^ and a* G Z~ , (97) 

where Ui^ are constants. The algorithm for the Painleve test is implemented by means of the following 
three steps: 

Step 1 (Determine the leading singularity or dominant behavior). We replace 

Ui (t) ~ u i<0 (t - t Q ) a * (98) 

into Q95JI to determine ct^ and it^o and we obtain an algebraic system with a^, assuming negative integer 
values and to arbitrary. 

We require that two or more terms of each equation may balance and determine a* and u^.o- 
If any ai is not integer, the system is not of Painleve type in its strong version. 

If there are more than one solution for ai or u^o they define branches and the following steps of the 
algorithm need to be applied for each of these branches. 
Step 2 (Determine the resonances). 

For each ai and u^.o, we calculate the integers r for which Ui iT is an arbitrary function in 1951 We 
replace the truncated series 

«i(t) = «i,o(* - to) a > + «i,r(* " to) a ' +r - (99) 
by (|95|l . and we look for integer r for which Ui_ r is an arbitrary constant . 

To do that, after replacing the truncated series by l|95|l . we keep the most singular terms in (t — to), 
and the coefficients of u iyr are set to zero. We get: 

Qu r = 0, Ur = (ui^ r U2. r ---UM,r) 1 \ (100) 

with Q an M x M matrix depending of r. 

The resonances are the roots of det(Q) = 0. 

In every system with the Painleve property, the resonance (—1) will be present and correspond to 
arbitrary (t — to). The resonance with zero value may also be present, depending of the number of 
arbitrary values u^o- 

Step 3 (Compatibility conditions and constants of motion). 

For every resonance found in the previous step, there is a compatibility condition which must be 
verified in order that the system pass the Painleve test. The compatibility conditions are verified by 
inserting 

ru 

u t (t) = (t-t ) a >J2 u *At-to) k (101) 

k=0 
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into (|95(l . where rjvf is the highest positive integer resonance. 

If all these compatibility conditions are satisfied so that they introduce a sufficient number of arbitrary 
constants, then the system is said to be of Painleve type. 
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